Extra-Chromosomal DNA Sequencing Reveals Episomal Prophages Capable of Impacting Virulence Factor Expression in Staphylococcus aureus

Staphylococcus aureus is a major human pathogen with well-characterized bacteriophage contributions to its virulence potential. Recently, we identified plasmidial and episomal prophages in S. aureus strains using an extra-chromosomal DNA (exDNA) isolation and sequencing approach, uncovering the plasmidial phage ϕBU01, which was found to encode important virulence determinants. Here, we expanded our extra-chromosomal sequencing of S. aureus, selecting 15 diverse clinical isolates with known chromosomal sequences for exDNA isolation and next-generation sequencing. We uncovered the presence of additional episomal prophages in 5 of 15 samples, but did not identify any plasmidial prophages. exDNA isolation was found to enrich for circular prophage elements, and qPCR characterization of the strains revealed that such prophage enrichment is detectable only in exDNA samples and would likely be missed in whole-genome DNA preparations (e.g., detection of episomal prophages did not correlate with higher prophage excision rates nor higher excised prophage copy numbers in qPCR experiments using whole-genome DNA). In S. aureus MSSA476, we found that enrichment and excision of the prophage ϕSa4ms into the cytoplasm was temporal and that episomal prophage localization did not appear to be a precursor to lytic cycle replication, suggesting ϕSa4ms excision into the cytoplasm may be part of a novel lysogenic switch. For example, we show that ϕSa4ms excision alters the promoter and transcription of htrA2, encoding a stress-response serine protease, and that alternative promotion of htrA2 confers increased heat-stress survival in S. aureus COL. Overall, exDNA isolation and focused sequencing may offer a more complete genomic picture for bacterial pathogens, offering insights into important chromosomal dynamics likely missed with whole-genome DNA-based approaches

Genetic evidence for the interaction between Bacillus anthracis-encoded phage receptors and their cognate phage-encoded receptor binding proteins

Bacteriophages such as γ and AP50c have been shown to infect strains of Bacillus anthracis with high specificity, and this feature has been exploited in the development of bacterial detection assays. To better understand the emergence of phage resistance, and thus the potential failure of such assays, it is important to identify the host and phage receptors necessary for attachment and entry. Using genetic approaches, the bacterial receptors of AP50c and γ have been identified as sap and GamR, respectively. A second AP50c-like phage, Wip1, also appears to use sap as a receptor. In parallel with this work, the cognate phage-encoded receptor binding proteins (RBPs) have also been identified (Gp14 for γ, P28 for AP50c, and P23 for Wip1); however, the strength of evidence supporting these protein–protein interactions varies, necessitating additional investigation. Here, we present genetic evidence further supporting the interaction between sap and the RBPs of AP50c and Wip1 using fluorescently tagged proteins and a panel of B. anthracis mutants. These results showed that the deletion of the sap gene, as well as the deletion of csaB, whose encoded protein anchors sap to the bacterial S-layer, resulted in the loss of RBP binding. Binding could then be rescued by expressing these genes in trans. We also found that the RBP of the γ-like prophage λBa03 relied on csaB activity for binding, possibly by a different mechanism. RBPλBa03 binding to B. anthracis cells was also unique in that it was not ablated by heat inactivation of vegetative cells, suggesting that its receptor is still functional following incubation at 98°C. These results extend our understanding of the diverse attachment and entry strategies used by B. anthracis phages, enabling future assay development

The META tool optimizes metagenomic analyses across sequencing platforms and classifiers

A major challenge in the field of metagenomics is the selection of the correct combination of sequencing platform and downstream metagenomic analysis algorithm, or “classifier”. Here, we present the Metagenomic Evaluation Tool Analyzer (META), which produces simulated data and facilitates platform and algorithm selection for any given metagenomic use case. META-generated in silico read data are modular, scalable, and reflect user-defined community profiles, while the downstream analysis is done using a variety of metagenomic classifiers. Reported results include information on resource utilization, time-to-answer, and performance. Real-world data can also be analyzed using selected classifiers and results benchmarked against simulations. To test the utility of the META software, simulated data was compared to real-world viral and bacterial metagenomic samples run on four different sequencers and analyzed using 12 metagenomic classifiers. Lastly, we introduce “META Score”: a unified, quantitative value which rates an analytic classifier’s ability to both identify and count taxa in a representative sample

A bioinformatic filter for improved base-call accuracy and polymorphism detection using the Affymetrix GeneChip® whole-genome resequencing platform

DNA resequencing arrays enable rapid acquisition of high-quality sequence data. This technology represents a promising platform for rapid high-resolution genotyping of microorganisms. Traditional array-based resequencing methods have relied on the use of specific PCR-amplified fragments from the query samples as hybridization targets. While this specificity in the target DNA population reduces the potential for artifacts caused by cross-hybridization, the subsampling of the query genome limits the sequence coverage that can be obtained and therefore reduces the technique's resolution as a genotyping method. We have developed and validated an Affymetrix Inc. GeneChip® array-based, whole-genome resequencing platform for Francisella tularensis, the causative agent of tularemia. A set of bioinformatic filters that targeted systematic base-calling errors caused by cross-hybridization between the whole-genome sample and the array probes and by deletions in the sample DNA relative to the chip reference sequence were developed. Our approach eliminated 91% of the false-positive single-nucleotide polymorphism calls identified in the SCHU S4 query sample, at the cost of 10.7% of the true positives, yielding a total base-calling accuracy of 99.992%

Genome Sequencing of Four Strains of \u3ci\u3eRickettsia prowazekii\u3c/i\u3e, the Causative Agent of Epidemic Typhus, Including One Flying Squirrel Isolate

Rickettsia prowazekii is a notable intracellular pathogen, the agent of epidemic typhus, and a potential biothreat agent. We present here whole-genome sequence data for four strains of R. prowazekii, including one from a flying squirrel

Transcriptional profiling with a blood pressure QTL interval-specific oligonucleotide array

Although the evidence for a genetic predisposition to human essential hypertension is compelling, the genetic control of blood pressure (BP) is poorly understood. The Dahl salt-sensitive (S) rat is a model for studying the genetic component of BP. Using this model, we previously reported the identification of 16 different genomic regions that contain one or more BP quantitative trait loci (QTLs). The proximal region of rat chromosome 1 contains multiple BP QTLs. Of these, we have localized the BP QTL1b region to a 13.5-cM (20.92 Mb) region. Interestingly, five additional independent studies in rats and four independent studies in humans have reported genetic linkage for BP control by regions homologous to QTL1b. To view the overall renal transcriptional topography of the positional candidate genes for this QTL, we sought a comparative gene expression profiling between a congenic strain containing QTL1b and control S rats by employing 1) a saturated QTL1b interval-specific oligonucleotide array and 2) a whole genome cDNA microarray representing 20,465 unique genes that are positioned outside the QTL. Results indicated that 17 of the 231 positional candidate genes for this QTL are differentially expressed between the two strains tested. Surprisingly, \u3e1,500 genes outside of QTL1b were differentially expressed between the two rat strains. Integrating the results from the two approaches revealed at least one complex network of transcriptional control initiated by the positional candidate Nr2f2. This network appears to account for the majority of gene expression differences occurring outside of the QTL interval. Further substitution mapping is currently underway to test the validity of each of these differentially expressed positional candidate genes. These results demonstrate the importance of using a saturated oligonucleotide array for identifying and prioritizing differentially expressed positional candidate genes of a BP QTL. Copyright © 2005 the American Physiological Society

Iterative microarray and RNA interference-based interrogation of the Src-induced invasive phenotype

Src kinase has long been recognized as a factor in the progression of colorectal cancer and seems to play a specific role in the development of the metastatic phenotype. In spite of numerous studies conducted to elucidate the exact role of Src in cancer progression, downstream targets of Src remain poorly understood. Gene expression profiling has permitted the identification of large sets of genes that may be functionally interrelated but it is often unclear as to which molecular pathways they belong. Here we have developed an iterative approach to experimentally reconstruct a network of gene activity regulated by Src and contributing to the invasive phenotype. Our strategy uses a combination of phenotypic anchoring of gene expression profiles and loss-of-function screening by way of RNA-mediated interference. Using a panel of human colon cancer cell lines exhibiting differential Src-specific activity and invasivity, we identify the first two levels of gene transcription responsible for the invasive phenotype, where first-tier genes are controlled by Src activity and the second-tier genes are under the influence of the first tier. Specifically, perturbation of first-tier gene activity by either pharmacologic inhibition of Src or RNA-mediated interference-directed knockdown leads to a loss of invasivity and decline of second-tier gene activity. The targeting of first-tier genes may be bypassed altogether because knockdown of second-tier genes led to a similar loss of invasive potential. In this manner, numerous members of a transcriptional cascade pathway for metastatic activity have been identified and functionally validated. ©2005 American Association for Cancer Research

Table_1_Extra-Chromosomal DNA Sequencing Reveals Episomal Prophages Capable of Impacting Virulence Factor Expression in Staphylococcus aureus.docx

<p>Staphylococcus aureus is a major human pathogen with well-characterized bacteriophage contributions to its virulence potential. Recently, we identified plasmidial and episomal prophages in S. aureus strains using an extra-chromosomal DNA (exDNA) isolation and sequencing approach, uncovering the plasmidial phage ϕBU01, which was found to encode important virulence determinants. Here, we expanded our extra-chromosomal sequencing of S. aureus, selecting 15 diverse clinical isolates with known chromosomal sequences for exDNA isolation and next-generation sequencing. We uncovered the presence of additional episomal prophages in 5 of 15 samples, but did not identify any plasmidial prophages. exDNA isolation was found to enrich for circular prophage elements, and qPCR characterization of the strains revealed that such prophage enrichment is detectable only in exDNA samples and would likely be missed in whole-genome DNA preparations (e.g., detection of episomal prophages did not correlate with higher prophage excision rates nor higher excised prophage copy numbers in qPCR experiments using whole-genome DNA). In S. aureus MSSA476, we found that enrichment and excision of the prophage ϕSa4ms into the cytoplasm was temporal and that episomal prophage localization did not appear to be a precursor to lytic cycle replication, suggesting ϕSa4ms excision into the cytoplasm may be part of a novel lysogenic switch. For example, we show that ϕSa4ms excision alters the promoter and transcription of htrA₂, encoding a stress-response serine protease, and that alternative promotion of htrA₂ confers increased heat-stress survival in S. aureus COL. Overall, exDNA isolation and focused sequencing may offer a more complete genomic picture for bacterial pathogens, offering insights into important chromosomal dynamics likely missed with whole-genome DNA-based approaches.</p